Higher-Order Exponential Integrators for Quasi-Linear Parabolic Problems. Part I: Stability

Explicit exponential integrators based on general linear methods are studied for the time discretization of quasi-linear parabolic initial-boundary value problems. Compared to other exponential integrators encountering rather severe order reductions, in general, the considered class of exponential general linear methods provides the possibility to construct schemes that retain higher-order accuracy in time when applied to quasi-linear parabolic problems. Employing an abstract framework, the considered problems take the form of initial value problems on Banach spaces: u (t) = Q(u(t))u(t), t ∈ (0, T), u(0) given. A fundamental requirement for the stability and error analysis is that the domains of the defining sectorial operators Q(v) : D = D(Q(v)) → X are independent of v ∈ V ⊂ X. The scope of applications in particular includes quasi-linear parabolic evolution equations subject to Dirichlet boundary conditions. The work is divided into two parts. In Part I, stability bounds in the norms of certain intermediate spaces between the domain D and the underlying Banach space X are deduced. In view of practical applications, the stability estimates are stated for variable time stepsizes, under mild restrictions on the ratios of subsequent stepsizes. The stability results provide a basic ingredient for the convergence analysis given in Part II. 1. Introduction. Our present work on efficient time integration methods for quasi-linear parabolic problems comprises two parts. In this first part, we introduce the considered class of variable stepsize explicit exponential general linear methods and study their stability behavior. A second part shall be concerned with the convergence analysis of these exponential integrators. Scope of applications. Quasi-linear parabolic initial-boundary value problems typically arise in the modeling of minimal surfaces and mean curvature flow, in the study of fluids in porous media and sharp fronts in polymers, and for the description of thin fluid films and diffusion processes with state-dependent diffusivity; see, for instance,, we may cast a quasi-linear parabolic